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ABSTRACT 

Context. The quiet Sun magnetic fields produce ubiquitous bright points (BPs) that cover a significant fraction of the solar surface. 
Their contribution to the total solar irradiance (TSI) is so-far unknown. 

Aims. To measure the center-to-limb variation (CLV) of the fraction of solar surface covered by quiet Sun magnetic bright points. The 
fraction is referred to as fraction of covered surface, or PCS. 

Mefhods. Counting of the area covered by BPs in G-band images obtained at various heliocentric angles with the 1-m Swedish Solar 
Telescope on La Palma. Through restoration, the images are close to the diffraction limit of the instrument (~ O.'l). 
Results. The PCS is largest at disk center 1 %), and then drops down to become =^ 0.2 % at /j ^ 0.3 (with /j the cosine of the 
heliocentric angle). The relationship has large scatter, which we evaluate comparing different subfields within our FOVs. We work out 
a toy-model to describe the observed CLV, which considers the BPs to be depressions in the mean solar photosphere characterized by a 
depth, a width, and a spread of inclinations. Although the model is poorly constrained by observations, it shows the BPs to be shallow 
structures (depth < width) with a large range of inclinations. We also estimate how different parts of the solar disk may contribute to 
TSI variations, finding that 90 % is contributed by BPs having /j > 0.5, and half of it is due to BPs with yu > 0.8. 
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1. Introduction 

Our understanding of the quiet Sun magnetic fields 
has drastically improved during the last decade 
(for recent reviews, see, e.g., de Wiin et al.1 l2009t 
ISanchez Almeida & Martinez GonzilezI 12011 ). We have 
gone from magn etic signals present on l y at the network 
boundaries (e.g., iBeckersI 119771 : ISolankil Il993h . to ubiqui- 
tous polarization sign als created thr ough Hanle effect (e.g., 
iFaurobert-Sc holl 1993; 'Tr uiillo Bueno et al. 2004) and Zeeman 

Sanchez Almeida & Lites 
2003 , ; ,Harvev et ah ,2007t 



(e.g., i Lin & Rimme ld Il999[ 



X 



■ effect 

' 12000*; 'Dommgue z Cerdena et al.l 
Lites et al. 2008). The wealth of quiet Sun magnetic structures 
makes them potentially import ant to understand the global 
magnetic properties of th e Sun ( ISanchez Almeida et alj 120031 
iTruiillo Bueno et al]|2004l) . and also makes it unhkely that the 
quiet Sun magne tism results from the decay of active regions 
(e.g., [S^nchez A lmeida et alj 12003). Theoretical arguments, 
corroborated by numeric al experiments, favor a different 
produ ction mecha nism (P etrovav & Szakalvl T993t cTttaneol 
19991: PVo gler & S chiissleJ l2007; "Pi etarila Graham et alj l2010t 
Moll et alj i201 li) . An efficient turbulent dynamo transforms 



into magnetic fields part of the kinetic energy of the granular 
convection. It generates a complex magnetic field which evolves 
in short time scales (a few min) and has small characteristic 
length-scales (< 1 Mm). 

In agreement with the turbulent dynamo scenario, quiet Sun 
magnetic fields come with strengths in the full range cover- 



ing from almost zero to 2 kG (Sanchez Almeida & Lites 200C; 
iDommguez Cerdeiia et al.''2b06; Martinez Gonzalez et al.,.2008 : 
Bommier et al. 2009; Viticchie et al. 201 1). Even if they only fill 
a small fraction of the quiet photosphere, the part having strong 
kG fields may be particularly important for a number of rea- 
sons. Firstly, the magnetic flux and energy increase with field 
strength, therefore, the energy and flux provided by kGs may 
surpass the contributio n of the more common but weaker fields 
dSanchez Almeida 2004). The need to co nsider kGs is illustrate d 
by the numerical experiments set up by ICameron et a 
where realistic granular convection redistribute initial hG fields 
so that daG, hG and kG field strengths have the same energy 
despite their very different area covering. Magnetic concentra- 
tions with kG fields may also be imp ortant because buoyancy 
makes them vertical (e.g., ISchussle31l986,) and so, they natu- 
rally provide a mechanical conne ction between the photosphere 
and the upper atmosphere (e.g., van Ballegooiien et al. 19981; 
Schriiver & Title 2003; Goodman 2004; Jendersie & Peter 200% 
fSanchez Almeida et al.l 120071) . They can function as guides 
that sustain magneto-acustic wave propagation, or be physical 
channels connecting plasmas of different atmospheric layers. 
Finally, kG magnetic concentrations are expected to be partic- 
ularly bright due to the so-called hot- wall effecn (.Spi'uit . 1 9771 : 



' The magnetic pressure suffices to maintain kG structures in me- 
chanical balance within the photosphere, therefore, they are evacuated 
and transparent, allowing us to look through into the sub-photosphere, 
which is generally hotter and so brighter. 
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ICarlsson et aP l2004t iKeller et alJ |2004 . They produce bright 
points (BPs) which, depending on the variation across the solar 
disk, and during the solar cycle, may even contribute to the Total 
Solar Irradiance ( TSI) variations as network and plage mag- 
netic fields do (e.g.. lLeanlll997l:lFrohlich & Leanl2004trFrohlichl 
l20Tll) . 

The finding of BPs in the quiet Sun was immediately iden- 
tified as the kG magnet ic concentrations infer red from po- 
larization measurements dSanchez Almeida et alJ 120 04;!). Their 
basic properties and their ubiquitous presence have been 
confirmed by a numbe r of re searchers (deWiin 
elet & Wiehr '2008'; 'sinchez Almeida 
_ 2009, 2O10j Goode et al. 201 
of minutes similar to that 



SI 



2008' 

Viticchie et al 



1 et alJ 
L et alT 



2005 



2007 



They vary on 

time scales of minutes similar to that of granulation, and 
many BPs are at the resolution limit of the current instrumen- 
tation (some 0.1 Mm). They are extremely common, filling 
at least 1% of the solar su rface and outnumbering granules 
dSanchez Almeida et al.ll20ra). 



Schnerr & Spruit! (1201 ih have recently estimated the excess 



of brightness produced by the quiet magnetic fields at the disk 
center, turning out to be of the order of 0.15 %. Their contribu- 
tion is larger than the 0.08 % TSI variations associated with cy- 
cle, however, determining the impact of these kG fields on TSI 
demands knowing their center-to-limb variation (CLV), as well 
as the variation of the quiet Sun magnetic fields with the solar cy- 
cle. These two properties, central to assess the role of quiet Sun 
fields on TSI, are poorly known. As far as the variation along the 
cycle is concerned, the claims in the literat ure are for little varia- 
tion if any ( Sanchez Almeida 2003; Shchu kina & Truiillo Buenol 
120031; iHarvevI 1201 0). Unfortunately, the unc ertainty of suc h 
claims can be as large as a factor two (Faurobert et al.ll200lh . 
There are several works on the CLV signals associated with the 
quiet Sun magnetic fields (e.g., iMartmez Gonzalez et a l] l2008t 
iLites et al. 2008; Lites 201 1!) however, to the best of our knowl- 
edge, nothing is known on the CLV of the quiet Sun magnetic 
BPs. This is precisely the subject of our work, i.e., providing a 
first observational description of the CLV of the quiet Sun BPs. 
To be more exact, we evaluate how the area covered by quiet Sun 
BPs varies with the position on the disk. 

Two main difficulties hinder the analysis. First, one has 
to use images with enough spatial resolution and of uniform 
quality, since the number of BPs depe nds critically on the 
spatial r esolution of the observations (iTitle & Bergej Il996t 
ISanchez A lmeida et al. 2010). This issue is sorte d out using data 
from the 1-m Swedish Solar Telescope (SST, IScharmer et alj 
12003 1) obtained in a single day during moments of excellent see- 
ing, and then restored to get images with uni form resolution 
close to the diffraction limit of the instrument ( van No ort et alj 
l2005h . Second, the BPs can b e misidentified with gra nule bor- 
ders and other structures (e.g., iBovelet & Wiehij[2007h . a prob- 
lem particularly severe in near-limb images. This second prob- 
lem is addressed resorting to the cumbersome method of eye- 
ball identification whi ch, however, has been p roven to be reli- 
able for our purpose (ISanchez Almeida et al.l [2010„ and refer- 
ences therein). 

The work is organized as follows; the observations are pre- 
sented in § |2l The actual CLV measurements are described in 
§ |3] In principle, such measurements has to be interpreted us- 
ing realistic models of magneto-co nv ection (such as tho se used 
for plages bv lCarlsson eTal]l2004l or lKeller et ai]|2004l) . where 
the 3-dimensional geometry of the photosphere is considered 
consistently. It would require repeating the analyses carried out 
on the observed images using synthetic images from the sim- 
ulations. Then the CLV coming from a comprehensive battery 



of numerical simulations should be compared with the observa- 
tions. This detailed realistic approach clearly exceeds the scope 
of the work. However, we attempt a toy-modeling which, despite 
its simplicity, considers the key ingredients that within the hot- 
wall paradigm determine the CLV (§|4li. The model is compared 
with the observed CLV in § 14.21 which constraints some of the 
properties of the magnetic structures. The implications for TSI 
variations are analyzed in § |5] These results are discussed and 
put into context in § |6] 



2. Observations and data reduction 

Quiet Sun inter-network regions away from active areas were 
observed on August 7, 2006, with the 1-m Swedish Solar 
Telesc ope (SST, Rogue de los Muchachos Observatory, La 
Palma; ls"charmer et al.l l2003', 2002). The data consist of time se- 
ries of images taken at eight different heliocentric angles, 6, with 
cos 6) ranging from 0.34 to 1. The logbook, in Table [T] in- 
cludes the heliocentric angle of each series as well as the UT of 
observation. Seeing was variable, with periods of excellence. As 
we explain below, we select for analysis only the best snapshots. 

Images were simultaneously recorded in three channels, two 
of them in G-continuum (A 4363.9 A; FWHM 11 A) yielding 
pairs of simultaneous in-focus and ou t-of-focus images for post- 
facto application of phase diversity (iPaxman et al.llI992l) . The 
third channel, in the G-band (A 4305.6A; FWHM 10.8 A), is 
the one used for our analysis. Each camera continuously gath- 
ered images at a rate of some 100 min"'. An additional fourth 
channel was devoted to Call H (A 3968.5 A; FWHM 1.1 A), 
which was used during observation to avoid active regions in the 
field of view (FOV). We employed two 10-bit Kodak Mega Plus 
1536 X 1024 cameras for G-continuum, a 12-bit Redlake Mega 
Plus II 1536 X 1024 camera for G-band, and a 10-bit Kodak 
Mega Plus 2024 x 2048 camera for Call H. The image scale 
was (y.'0405pix"' in all cases. 

After dark-current subtraction, flatfielding, and elimination 
of spurious pixels and borders, the image restoration in the 
first three channels was performed simultaneously employ- 
ing the Multi -Object Multi-Frame B lind Deconvolution method 
(MOMFBD; Ivan Noort et"an l2005h . Sets of -30 images per 
channel (i.e., a total of ~ 3 x 30 images) were combined to pro- 
duce a single pair of simultaneous G-band and G-continuum re- 
stored snapshots. Thus, we achieved two time series of restored 
images with a cadence of ~10 s, an effective FOV of 58'.'6 x 
38'.'6, and an angular resolution close to the diffraction limit of 
the SST at the working wavelengths (~0'.'l). The images of these 
series were de-rotated to compensate diurnal field rotation, rigid- 
aligned, destretched to remove image distortion, and subsonic 
filtered to suppress p-modes and residual jitter stemming from 
destretching (Title et al. 1989). 

The final products of the reduction process were two time 
series (G-band and G-continuum) for each one of the eight he- 
liocentric positions on the solar disk, with a duration ranging 
from 28 to 100 images (from 5 to 17 min) as summarized in 
Table [T] They represent a set having homogeneous angular res- 
olution and covering a wide range of heliocentric angles. Based 
on their contrast, we select the G-band image of best quality 
around the middle of each time series. They are the reference 
images used in our CLV analysis. Examples of two such images 
are shown in Fig.[T] 
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Table 1. Description of the data sets including the fraction of covered surface FCS. 



Series 


Initial time 


Final time 


Duration 


# of images 








C" 


1 

1 


o. 1 






1 no 

iUU 


n QQQ -t- 


n no? 
yj.yjyjz. 




1 ozi -t- n 1 ^ 


2 


8:47:16 


8:54:36 


0:07:20 


42 


0.610 ± 


0.028 


0.36 + 0.13 


1.16 + 0.18 


3 


8:56:41 


9:12:58 


0:16:17 


98 


0.582 ± 


0.030 


0.28 + 0.13 


1.21 +0.15 


4 


9:16:12 


9:29:02 


0:12:50 


78 


0.341 + 


0.065 


0.19 + 0.07 


1.15+0.13 


5 


9:32:44 


9:48:54 


0:16:10 


96 


0.802 ± 


0.027 


0.83 ± 0.20 


1.10 + 0.15 


6 


9:49:53 


10:06:25 


0:16:32 


98 


0.995 ± 


0.004 


0.88 + 0.10 


1.05+0.15 


7 


10:09:02 


10:23:13 


0:14:11 


86 


0.926 ± 


0.013 


0.58 + 0.15 


1.13+0.16 


8 


10:27:03 


10:31:50 


0:04:47 


28 


0.527 ± 


0.043 


0.31+0.10 


1.22 + 0.15 



" Mean intensity and rms variation relative to mean photosphere 



3. Center-to-limb variation of the solar surface 
occupied by bright points 

We have measured the fraction of solar surface occupied by BPs 
in each one of the eight G-band reference images described in 
§ |2] which correspond to different heliocentric positions on the 
disk. From now on this quantity will be referred to as Fraction 
of Covered Surface (FCS). 

The detection of BPs has been performed by eye with the 
help of com puter tool s. The procedure is the same used by 
[Sanchez Almeida et alj ( |2004 l20ldl) . and can be briefly de- 
scribed as follows: the reference im age of a part icular time series 
is segmented with the algorithm bv lStrousI (1199 4) that, based of 
the sign of the second derivative of the intensity at every im- 
age point, detects areas locally bright (i.e. areas brighter than 
their surroundings). The algorithm creates a binary mask with 
the pixels of the bright patches set to 1 . Using an interactive pro- 
gram, we flick on the computer display the reference image and 
the binary mask which allows us to identify coincidences be- 
tween G-band BPs and small segmented patches. The BPs are 
selected one by one, considering as such those matching a seg- 
mented patch in the binary mask when they overlay intergranular 
lanes and preserve identity along several frames close to the ref- 
erence image (at least plus minus two frames). We also select 
some faint bright features that are not evident BPs in the refer- 
ence image but which reveal themselves as such in the preceding 
or the following image. As a general policy doubtful structures 
are discarded. The procedure is carried out for two different sat- 
uration levels of the reference image on the computer display - 
low contrast to select the brightest and more evident BPs, and 
high contrast to identify weaker structures. 

The area occupied by the BPs is calculated as the number 
of pixels corresponding to the selected BPs. The ratio between 
this area and the total number of pixels in the image provides 
the FCS. Thus, we get the FCS for every yU, as listed in Table [1] 
These values should be regarded as lower limits because of three 
reasons: (1) we only include secure BPs, (2) the number of detec- 
tions increases with the angular resolution, which is finite, and 
(3) the segmenta tion algorithm underestimates the area covered 
by the large BPs dSanchez Almeida et al.ll20l"oh . As an argument 
for consistency, we point out that the FCSs at disk center are in 
agreement with that obtained from images of similar quality us- 
ing the same method bv lSanchez Almeida et alJ (1201 Ol ~ 0.9%), 
but are higher than those obt ained from lower quality images 
dSanchez Akneida et allllool -0.5%). 

The CLV of the FCS values are displayed in Fig. |2] Each 
symbol corresponds to one of the series, and includes error bars 
both in heliocentric angle and FCS. The horizontal error bars 
represent the range of // in the FOVs. The vertical error bars have 



been estimated dividing each full FOV into 24 non-overlapping 
subfields. They represent the standard deviation of the mean 
FCS, considering the various FCSs in the independent sub-fields. 
These error bars estimate the statistical error of the measure- 
ment. They do not account for the systematic errors stemming 
from the subjectivity of our BP identification. 

4. Interpretation 

4.1. A toy-model kG magnetic concentration 

In essence, a kG magnetic concentration represents a depression 
in the solar surface. Light comes from deeper photospheric lay- 
ers which are usually hotter, and so the structure looks brighter. 
This basic idea, simplified to the extreme but retaining key in- 
gredients, is considered in the 2-dimensional toy-model of mag- 
netic concentration represented in Figure [3^. It portrays the ver- 
tical section of a magnetic concentration of width a, which is 
located at an heliocentric angle 6. Due to the evacuation, we see 
deeper through the magnetic concentration, and the difference 
of geometrical depths with respect to the non-magnetic photo- 
sphere is denoted by h. The concentration can be inclined with 
respect to the local vertical, an inclination that we parameterize 
as ip. We assume the magnetic concentration to look bright only 
if our line-of-sight (LOS) reaches its bottom. In other words, the 
fraction of solar atmosphere that looks bright according to this 
model is proportional to I (see Fig. [3^). For a given configuration 
of the concentration (set by a, /i, and tp), it shows up bright (/ 9^ 0) 
only for a limited range of heliocentric angles, i.e., & < 6 < & 
as illustrated in Fig. [3j5. These extreme heliocentric angles are 
given by, 

a a 

tan0 = tan^H-- ; tan© =tani^-- , (1) 

h cos (f h cos 

with the natural constrains < © < 90° and < ©' < 0. Given 
an heliocentric angle within the proper interval (0' ,©), the 
length of the concentration that is visible, I, turns out to be, 

/ = — /z I tan - tan I . (2) 

cos if 

The FOV contains a number of these idealized magnetic con- 
centrations having different properties, therefore, within this toy 
model, the FCS is proportional to the mean / averaged over the 
set of magnetic concentrations, {I), 

FCS =/(/). (3) 

The scaling factor / accounts for the number of concentrations 
per unit surface, the area of the individual concentrations, and 
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5 10 lb 20 2b 




b 10 lb 20 2b 



Fig. 1. Images to illustrate how the quiet Sun G-band BPs look 
at different heliocentric angles - top ji - 0.99, bottom ji - 0.53. 
They correspond to the reference images in the series number 
6 and 8, respectively (Table [TJ. Axes represent arcsec from the 
lower left corner. 



possibly other factors of order unity associated with the fact that 
the magnetic concentrations do have a 3 -dimensional structure 
not included in our 2-dimensional toy-model. In our case we ad- 
ditionally assume all concentrations to have the same a and h, 
but different inclination^ with the tilt ip following a uniform 
distribution from -ipo to +ip() ((/jq < 90°). Consequently, our 
model for the CLV of FCS depends of three parameters, namely, 

^ We consider a distribution of tilts because the simpler alternative of 
assuming purely vertical magnetic concentrations produces a CLV with 
a sharp drop at small /j, which is not observed. 




0.0 0.2 0.4 0.6 0.8 1.0 
= cos d 



Fig. 2. Measurements of the fraction of solar surface occupied 
by BPs, FCS, at 8 heliocentric positions. Vertical bars represent 
the measure errors while horizontal bars show the range of he- 
liocentric positions included in the finite extension of the FOV. 
The three different curves represent the three different solutions 
provided by our toy-model - the solid, the dashed and the dot- 
dashed lines correspond to models a, b, and c in Table|2] respec- 
tively. 




(a) (b) 

Fig. 3. (a) Vertical section of the toy-model magnetic concen- 
tration used in the paper (the thick solid line). It represents a 
concentration of width a, tilted an angle ip with respect to the 
local vertical and observed at heliocentric angle Q. The symbols 
h and / are, respectively, the geometrical depth of the concentra- 
tion and the portion of the bottom accessible to observation, (b) 
For a given tilt angle ip, and 0' represent the extreme inclina- 
tion angles that allow us to observe the bottom of the magnetic 
concentration. The encircled plus symbol between figures shows 
the positive sense of rotation used in our equations. 



hi a, ipQ and af. It does not depend on h and a separately since 
the scaling factor af absorbs the dependence on a of Z (see equa- 
tion 131). The closer to th e limb, the higher in the photosphere we 
observe (e.g. Stixlll99ll) . a fact that does not contradict our as- 
sumption of /j/a to be independent of heliocentric angle. We are 
assuming that this systematic rising towards the limb affects the 
photosphere globally. In other words, both magnetic concentra- 
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tion and surroundings are uplifted by the same amount, leaving 
unchanged the relative depth of the concentration h. 

Obviously, it is difficult to justify the use of our toy-model 
based on its realism. The Sun is not 2-dimensional, the magnetic 
concentrations are not slabs, their walls are also bright, the non- 
magnetic background is not uniform, several concentrations may 
contribute along a single LOS, the depth h depends on the field 
strength, and so on and so forth. However, the model consid- 
ers the basic physical ingredients responsible for the CLV of the 
FCS, i.e., the bright structures are depressed with respect to the 
mean photosphere, so that they can be observed or not depending 
on their width, depth, and the relative orientation between their 
axes and the LOS. Moreover, the use of the proper 3-dimentional 
numerical models of magneto-convection is complicated, and 
clearly goes beyond the scope of the work (see §[T]i. Using mod- 
els that consider 3-dimensional magnetic flux-tubes fanning- 
out with height may be viewed as an intermediate alternative, 
but we disregarded this possibility from scratch since the ana- 
lytic model become extremely involved, and its realism remains 
way off that of the proper numerical simulations of magneto- 
convection (cf. Buente et al. 1993| with |Carlsson et al. 2004 or 
IkeTler et al.ll2004ft . The use of simplified models is also justified 
by the large scatter in the observed CLV (Fig.|2]i, provided that 
they only intend to explain general trends. 



A 
V 



1.0 - 



0.8 



0.6 



0.4 - 



0.2 



0.0 - 




0.0 



0.2 



0.4 



0.6 
= cos 6 



0.8 



1.0 



Fig. 4. Fraction of magnetized surface producing BPs. The three 
curves coiTespond to the predictions of our toy-model for the 
three solutions portrayed in Fig.|2] Note that magnetic structures 
no longer render BPs as one approaches the limb. 



4.2. Results of the model 

As we pointed out in the previous section, the CLV in our toy- 
model depends on only three independent parameters, that we 
choose to be h/a,(p() and af. In order to determine which set 
of parameters reproduces the observation best, we caiTy out a 
least squares fit to determine the three free parameters. Since the 
model is non-linear, the fit has to be carried out using an iter- 
ative procedure for which we^mploy the standard Levenberg- 
Marquardt algorithm (e.g.. press et al.lll988i) . Non-linear least 
squares fits do not grant uniqueness. The algorithm seeks and 
finds local minima of the merit function (x^), and if it has sev- 
eral, the best fitting parameters depend on the starting point. In 
order to get rid of this undesired dependence, we repeated the 
fit using 10** different initializations, where the starting h/a,(pi) 
and a f were obtained from three independent uniform distribu- 
tions spanning the full range of sensible values (0.2 < h/a < 4, 
0° < (fo < 89° and 5 x lO^^ < af < 5 x lO"'). The ensemble 
average needed to account the distribution of tilts from -(po to (po 
(see § 14.11 ) was evaluated numerically from lO'* samples drawn 
from a uniform distribution. Essentially, the different initializa- 
tions converge to three different solutions. Table |2] summarizes 
the parameters characteristic of these three families - it gives 
the mean values among all solutions in the family plus the for- 
mal error bars provided by the -minimization algorithm. The 
corresponding fits are plotted in Fig. [2] The most recurrent result 
(73 % of the accepted solutions - a in Table |2]l is also the most 
reliable since it represents the smallest close to 1 (Table |2]i. 
It corresponds to fairly shallow magnetic concentrations (aspect 
ratio h/a ~ 3/5) with a large range of tilts (ipo w 64°). The sec- 
ond most common solution (b in Table |2]i also corresponds to 
shallow magnetic concentrations (h/a 3/10) but they are ver- 
tical. Its is slightly larger than that for solution a. The third 
case is intermediate between a and b. 

According to the model and its assumptions, one just sees 
as BPs a portion I of the surface, which is only a fraction of the 
true surface occupied by magnetic structures A - see Fig. [3] The 
two quantities agree when the LOS coincides with the axis of 
the concentration, but in general A > /. In our toy-model the 



fraction of true surface that is observed as BP turns out to be 
(/)/(A), where the brackets account for the fact that we consider 
an ensemble of magnetic concentrations. The variation with he- 
liocentric angle of this fraction is represented in Fig.|4] As one 
approaches the limb, most of the magnetic concentrations are no 
longer BPs - less than 20 % of them at yu < 0.5. Figure |4] also 
shows how the solution with large range of magnetic field incli- 
nations (a in Table |2] and the solid line in Fig.|2| hides a large 
fraction of the existing magnetic structures; only 60% of them 
are observed as BPs at disk center 

The closer to the limb the more difficult the BP identification 
is (§ |2]and Fig. [TJ. In order to explore the effect of this uncer- 
tainty, we repeated the fits removing the two FOVs closest to the 
limb (/i = 0.34 and 0.53). The results are shown in Table[3]and 
Fig- El Table[3]reflects the presence of four families of solutions, 
all of them with similar The third most frequent result (f in 
Table[3]- 17 % of cases) depicts a shallow flux tube h/a x 3/5 
with a significant tilting range (po ^ 46°, i.e., it looks quite simi- 
lar to the main solution a in Table|2] However, the most common 
solution (40 % of cases - d in TableO represents a much deeper 
tube, h/a > 3, although with a similar tilting range (i^o ~ 52°). 
Finally, the solution e in Table [3] resembles very much to that 
reported as solution c in Table |2l Removing the two most un- 
certain FCSs yields solutions similar to the full case, but it also 
introduces another possibility of rather slim magnetic concentra- 
tions (d in Table [3]). 

In short, the large scatter of the observations makes it easy 
to reproduce the CLV with our toy-model. The most favored so- 
lution coiTesponds to BPs being rather shallow structures with a 
large spread of vertical inclinations. This solution is not unique, 
and the model also allow for predominantly vertical structures 
and for narrow concentrations. 

5. Fraction of solar disk covered by BPs 

FCS(ju) describes the fraction of surface covered by BPs at each 
location on the Sun. Once FCS(/i) is known, the fraction of the 
solar disk covered by BPs can be determined by integration. The 
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Table 2. Parameters of the best fitting toy-model considering the eight heliocentric angles 



X' hi a ifio [°] af < A> solutions [%] 

0.90-0.91 0.65 ±0.16 64 + 5 (9.7 ± 1.6) x 10^^ 65.9 a 72.8 

1.45-2.00 0.29 ±0.01 1 ± 44 (8.0 ± 0.3) x 10"' 50.0 b 23.5 

1.97-1.98 0.43 ±0.02 6 ± 24 (8.8 ± 0.3) x 10-^ 50.1 c 2.8 



Table 3. Parameters of the best fitting toy-model considering the six center-most fields 





hi a 






af 


< A > 


solutions 


[%] 


0.47 


3.3 ± 1.6 


52 ± 1 


(4.7 


± 2.4) X 10-2 


59.0 


d 


39.5 


0.51-0.52 


0.46 ± 0.04 


7 ±29 


(8.9 


± 0.4) X 10-3 


50.1 


e 


29.9 


0.54 


0.59 ± 0.24 


46 ±33 


(9.7 


± 1.6) X 10-3 


56.6 


f 


17.2 


0.52 


1.11 ±0.36 


54 ±5 


(1.5 


± 0.5) X 10-2 


59.8 


g 


10.9 



Co 



0.010 



0.008 



0.006 



0.004 



0.002 



0.000 - 



0.0 




0.4 0.6 
/u. = cos 6 



Fig. 5. Same as Fig. |2]but for fits to the six center-most data 
points describing the CLV. The solid line, the dashed line, the 
dot dashed line, and the triple-dot dashed line correspond to so- 
lutions d, e, f and g in Table [3] respectively. 



fraction of solar disk with heliocent ric angle between /i and /i 
dji is Ifidfi (e.g.. lFoukal et al.lll99lh . therefore. 



Jo 



(4) 



provides the fraction of surface with BPs from the limb to fi. 
Obviously, F{1) represents the fraction covering the full Sun. 
F(fi) is given in Fig. |6j3 for the three solutions in Fig. |2l The 
three functions are very similar with F{1) ^ 0.005, i.e., some 
0.5 % of the solar disk is covered by BPs. The disk center is the 
maximum contributor to this BP covering - 90 % of it is due 
BPs at yu > 0.5, and half of the covering is associated with BPs 
at yu > 0.8 (Fig. I6J3). Even though the above conclusions use 
the CLV of our toy-model, it is important to emphasize that they 
are independent of the assumptions we made to work it out. The 
model CLVs are used here as smooth continuous approximations 
to describe the actual observations. 

As we argue in the previous paragraph, most of the area oc- 
cupied by BPs is at the disk center It seems that any potential 
contribution of the quiet Sun magnetic fields to the TSI cannot 
come from the solar limb. In order to make this statement more 



quantitative, we follow the equations bv'Fouka l et alj (1199 Ih to 
write down the contribution to the TSI of the BPs, i.e., it is given 

as 

with 

H(m) FCS(ix') [C(ti') - l] (3/z' + 2)fi' dfi'. (5) 

The symbol AS stands for the change on the solar irradiance S 
produced by the quiet Sun BPs. The previous equation assumes 
an Eddington limb-darkening for the quiet Sun intensity. The 
parameter C(jf) is the contrast of the BPs relative to the mean in- 
tensity averaged over all wavelengths. (Just to make the meaning 
of AS IS more clear: if BPs twice as bright as the mean photo- 
sphere fully cover the disk, i.e. C = 2 V yu, then AS IS = 1.) 
As we explain above, the part of the integrand of equation Q 
involving FCS is strongly biased towards the disk center. The 
only parameter that may counter-balance this effect is C. If C 
increases too rapidly with fi it may dominate the integrand, and 
so, the TSI. We ignore how C varies from center to limb, how- 
ever, as an educated-guess to evaluate H{jj), we have assumed 
C to scale as the intensity variation we measure in the G-band. 
This CLV of the contrast in the G-band is given in Table [T] and 
corresponds to the symbols represented in Fig.|6t. Two smooth 
curves fitted to the observed points are also included in the fig- 
ure - they consider the two extreme scenarios of the contrast re- 
maining constant at the limb (the dotted line) or dropping down 
significantly (the triple-dot dashed line). 

Figure|6}; displays H(fj.) for the three solutions in Fig.|2]con- 
sidering the contrast decreasing at the limb. They are very sim- 
ilar showing in all cases that the irradiance is strongly biased 
towards the disk center. This result is independent on whether 
the CLV of C has or not the drop at the limb - the thick solid 
line in Fig.|6]; is equivalent to the thin solid line except that the 
used C maintains a constant value at the limb (the dotted line in 
Fig. |6ji). From H(fi) in Fig. |6]: one finds that 90 % of the TSI is 
contributed by BPs having /i > 0.5, and some half of it is due to 
BPs with ^ > 0.8. 

The peak value of H is about 0.07 % (Fig.|6};). It corresponds 
to BPs with a typical contrast of the order of 1.15 filling the 
observed disk coverage of some 0.5 % (Fig.|6]3). This figure for 
the contribution of the qu iet Sun magnetic fields to TSI is similar 
to the recent estimate by ISchnerr & Spruill (1201 Ih mentioned in 
§[T] Actually, it is a factor of two smaller, but their estimate refers 
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0.0004 



0.0002 



0.0000 



Fig. 6. (a) Contrast of the BPs in the G-band as observed in our 
fields, C. The symbols are observations, and the lines represent 
two different smooth fits to the observed CLV. (b) Fraction of 
total solar disk covered by BPs from the solar limb to a particular 
fi-F in equation (|4]i. The three different types of line correspond 
to the three curves in Fig. |2] (c) Fraction of TSI contributed by 
the quiet Sun BPs - H in equation (|5]l. The three different types 
of thin lines correspond to the three curves in Fig. |2] with the 
contrast given by the triple-dot dashed line in (a). The thick solid 
line is equivalent to the thin solid line except that the contrast is 
taken as the dotted line in (a). 



to the disk center observed at a different wavelength. In view of 
the uncertainties involved in this type of work, the coincidence 
is worth pointing out. The two estimates are completely different 
and yet provide consistent results. 



6. Discussion and Conclusions 

We have measured the center-to-limb variation (CLV) of the area 
covered by G-band bright points (BPs) in the quiet Sun (fraction 
of covered surface or FCS). It is a parameter difficult to deter- 
mine since the detection of BPs critically depends on the angular 
resolution of the observation (§[Tli. We employ several time se- 
ries taken in two hours during moments of excellent seeing with 
the SST. The images were post processed using MOMFBD (see 
§ |2]l which provides a homogeneous set of images adequate for 
these subtle measurements. They were restored to provide an an- 
gular resolution close to the diffraction limit of the instrument at 
the working wavelength (some (X'l at the G-band). 

We find the FCS to be largest at disk center (^ 1 %), and 
then it drops down to become ^ 0.2 % at yu ^ 0.3. The relation- 
ship has large scatter, which we managed to estimate comparing 
different subfields within our FOVs (see the error bars in Fig. |2] 
and Table [U. The value obtained at the disk center agrees with 
previous estimates based on data of similar quality. 

We work out a toy model to describe the observed CLV. It as- 
sumes the magnetic bright points to be depressions in the mean 
solar photosphere, characterized by a depth, a width, and with a 
spread of inclinations. It is only an exploratory modeling which, 
however, includes the physical ingredients that seems to be re- 
sponsible for a kG magnetic concentration to show up as a bright 
feature on the solar disk. The solutions offered by our toy-model 
are poorly constrained, but they seem to show the BPs to be shal- 
low structure^ (ratio depth to width ^ 0.7 + 0.2) with a large 
spread of inclinations (^ +70°). 

Among others, the FCS is of interest because it determines 
the impact that quiet Sun magnetic fields may have on TSI vari- 
ations, an influence so far unknown. Since the measured FCS is 
so peaked towards disk center, any role that quiet Sun magnetic 
fields may have to play on TSI will be due to the center of the 
disk. According to our estimate, 90 % of the TSI is contributed 
by BPs having /i > 0.5, and half of it is due to BPs with /i > 0.8 
(§ EJl. This estimate is based on assuming the CLV of the BP 
(wavelength integrated) flux to scale as the observed G-band in- 
tensity. It is an ad-hoc assumption adopted for lacking of a better 
one which, however, has allowed us to provide a first constrain 
on the effect of quiet Sun BPs on TSI. 

In this sense, we have to stress that the FCS has been mea- 
sured in the G-band, because the magnetic concentrations are 
particularly conspicuous at this wavelengths. Then we implicitly 
assume throughout the work that the FCS is the same at all wave- 
lengths, i.e., that the area covered by magnetic concentration 
does not change with wavelength. (The G-band BPs are brighter 
but not bigger) This assumption remains to be proven however, 
for the time being, it seems to be a reasonable working hypothe- 
sis. 

Our toy-model suggests the BPs to be shallow structures 
with varied inclinations. One could test these predictions ob- 
serving polarization signals of quiet Sun BPs at disk center. On 
the one hand, the polarization signals provide magnet ic field in- 
clinations through Stokes inversion (e.g., ISocas-Navarro et al.l 
12008 ). On the other hand, assuming the magnetic concentra- 



' This result is not inconsistent with the magnetic concentra- 
tions producing BPs being modeled as a compactly-packed ensem- 
ble of narrow magnetic concentrations arranged in a micro- structured 
magnetic atmosphere MI SMA fashion t Sanchez. Almeida et al.lll993 : 
Sanchez Almeida 2000; S anchez Almeida et alj|2001l) . The CLV of the 
FCS is sensitive to the global scale of the ensemble, whereas the 
MISMA accounts for the smallest scales responsible, among others, for 
the asymmetries of the spectral lines formed in the atmosphere. 
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tions to be in mechanical balance, o ne can infer their depths 
(e.g. JSanchez Almeida & Litesll2000h . These inferences require 
low-noise spectro-polarimetry with an angular resolution sim- 
ilar to that of our G-band images. They represent a techni- 
cal challenge, but suc h observations se em to be doable in the 
near future (see, e.g., iLagg et al]|2010l) . As we mentioned in 
the introduction, the existing numerical simulations of magneto- 
convection predict the BPs to b e depressed wit h respect to the 
mean photosphere (e.g.. lVogler e t al. 2005; Vo gler & Schiissle^ 
L2007"). Whether the faint BPs observed in quiet Sun are predicted 
to be superficial or deep remains to be worked out. However, the 
existing simulations of plage magnetic concentrations suggest 
the continuum intensity to be formed in a rather shallow region 
(iKeller et al.ll2004tlCarlsson et al.ll200l . 
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